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P-mode observations on ex Cen A^ 
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Abstract. We have made a clear detection of p-mode oscillations in the nearest solar-like G2V star a Cen A 
with the CORALIE spectrograph on the 1.2-m Swiss telescope at the ESO La Silla Observatory. We report the 
5 nights of observation on this star in May 2001 during which 1260 high precision radial velocity measurements 
were obtained. The power spectrum clearly shows several identifiable peaks between 1.7 and 3 mHz. The average 
large splitting of 105.7 fiHz and the amplitude of about 35 cm s~^ of these modes are in agreement with theoretical 
expectations. 
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1. Introduction 

One of the major objectives of the stellar physic is the 
understanding of stellar interiors and the modelisation of 
star evolution. Asteroseismology, which consists to mea- 
sure amplitudes and frequencies of oscillation modes is 
an ideal tool to constrain stellar models and evolutionary 
theory. 

Recent results of ground-based observations using 
Doppler techniques have permitted to detect solar-like os- 



cillations in the F5IV-V star Procyon (Martic et al. 1999) 
andthe G2IV star /? Hydri (Bedding et al. ^00l| Carrier et 
al. 2001). These two stars exhibit good evidence for excess 
power centered at 1 mHz with peak amplitude of about 
50 cms~^. 

A primary and challenging target for the search for 
p-mode oscillations is the G2V star Rigil Centaurus 
{a Cen A, HR5459). Thanks to its proximity and its 
belonging to a binary system, its characteristics are 
well-determined which simplifies the interpretation of 
asteroseimological results. Scaling from the solar case, 
the frequency of its greatest mode is expected to be 
i^max = 2.3 mHz, the primary frequency splitting Ai'o = 
105.8 //Hz and the oscillat ion am plitude Aqsc = 31.1 
cms~^ (Kjeldsen & Bedding 1995a ). Detailed Models with 
seismological analysis have been rec ently proposed for 
a Cen A (see, e.g., Pourba ix et al. 1999 , Guenther & 
Demarque |2000| , Morel et al. |2000D . 
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Based on observations collected with the CORALIE echelle 
spectrograph on the 1.2-m Euler Swiss telescope at La Silla 
Observatory, ESO Chile 



Several groups have already conducted thorough at- 
tempts to detect signature of the p-mode oscillations on 
this star. Two groups claimed mode detections with am- 
plitudes 3.2-6. 4 gre ater than solar (Gelly et al. 1986 ; 
Pottasch et al. 1992). It was however infirmed by three 



others groups which obtained upper limits of mode am- 

solar (Brown I 
Kjeldsen et al 



Edmonds & Cram |1995|; Kje ldsen et al. |1999|). More re- 
cently Schou & Buzasi (2001) have made photometric ob- 



plitudes of 1.4-3 t imes solar (Brown & GiUiland 1990 
|1995|; 



servations of a Cen A with the WIRE spacecraft and re- 
ported a possible detection of P-modes. 

We obtained high precision radial velocity (RV) mea- 
surements on a Cen A and report in this paper our obser- 
vations and the power spectrum analysis of this star. 

2. Observations and data reduction 

The observations of a Cen A were carried out with the 
CORALIE fiber-fed echelle spectrograph (Queloz et al. 
^000 ) mounted on the 1.2-m Swiss telescope at the ESO 
La Silla Observatory. CORALIE is the southern hemi- 
sphere twin of the ELODIE spectrograph (Baranne et 
al. 1996), both of them well known for their discover- 



ies of extra-solar planets. The wavelength domain ranges 
from 3875 to 6820 A recorded on 68 orders. Thanks to 
a slightly different optical combination at the entrance of 
the spectrograph and the use of a 2k by 2k CCD with 
15-fim pixels, CORALIE reaches a spectral resolution of 
50 000 with a 3 pixels sampling. The total efficiency in- 
cluding atmosphere, telescope, spectrograph and detector 
was measured and reaches about 1.5 % at 5500 A in the 
best seeing case. During stellar exposures, the spectrum of 
a thorium lamp carried by a second fiber is simultaneously 
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Table 1. Distribution and dispersion of RV measurements 






••y. 
w 










v^. 


• * • 


i 

»v i 

• *1 






"•J I* 

If 



JD - 2452029 



Fig. 1. Radial velocity measurements of a Cen A. A 3 or- 
der polynomial fit subtraction was applied on each night. 
The dispersion reaches 1.47 ms^^. 



recorded in order to monitor the spectrograph instability 
and thus to obtain high precision velocity measurements. 
Spectra are extracted and reduced in real time using the 
INTER- TACOS software package developed by D. Queloz 
and L. Weber at the Geneva Observatory (see Baranne et 
al. 



1996) 



The standard RV computation by a cross-correlation 
algorithm was replaced by an algorithm based on the op- 
timum weight procedure first proposed by Connes ( 1985| ). 
We developed this method and showed that the optimum 
procedure is a factor 1.6 more efhcient in term of pho- 
ton noise than the cross-correlation procedure (Bouchy et 
al. ^001 ). The method uses the full spectral information 
available and compares two spectra of the same source to 
compute their velocity difference. Its present limitation is 
that the Doppler shift must remain small compared to the 
line-width. This is the case in the time scale of few hours 
when the velocity change, dominated by Earth motion, 
is less than 100 ms~^ (~ 1/50 of a typical stellar line- 
width). During a whole night, Earth motion can introduce 
a Doppler shift greater than 500 ms~^. An error of about 
3 to 5 % is then introduced if an unique "reference" spec- 
trum is used during the whole night. However, this effect 
can be entirely eliminated by a low order polynomial fit 
subtraction without consequences for the characterization 
of high frequency oscillation modes. 

a Cen A was observed over 5 nights in May 2001. A 
journal of these observations is given in Table 1. We took 
sequences of 40-s exposures with a dead time of 110-s in- 
between. In total, 1260 spectra were collected with typical 
signal-to-noise ratio in the range 300-420 at 550 nm. 



Date 


Nb spectra 


(J ( ms ^) 


2001/04/29 


220 


1.58 


2001/04/30 


259 


1.49 


2001/05/01 


268 


1.39 


2001/05/02 


246 


1.40 


2001/05/03 


267 


1.53 



Radial velocities are computed for each night relatively 
to the highest signal-to-noise ratio spectrum obtained in 
the middle of the night. The resulting velocities are pre- 
sented in Fig. 1. The dispersion of these measurements 
reaches 1.47 ms~^ and the individual value of each night 
is listed in Table 1. The precision of our RV measure- 
ments is compared with the fundamental uncertainty due 
to photon noise using the method described by Bouchy 
et al. (2001). The uncertainties coming from the stellar 
spectrum and from the simultaneous thorium spectrum 
used in the instrumental tracking are typically the same 
and equal to 0.50 ms""'^. The quadratic sum of these two 
photon noise contributions is 0.70 ms^^. 

3. Power spectrum analysis 

In order to compute the power spectrum of the velocity 
time series of Fig. 1, we used the Lomb-Scargle modified 
algorithm (Lomb 1976 , Scargle 1982 ) for unevenly spaced 
data. The resulting LS periodogram, shown in Fig. 2, ex- 
hibits a series of peaks between 1.7 and 3 niHz modu- 
lated by a broad envelope, which is the typical signature 
of solar-like oscillations. This signature also appears in 
the power spectrum of every individual night. Toward the 
lowest frequencies {v < 0.5 mHz), the power rises and 
scales inversely with frequency squared as expected for 
instrumental instabilities. The mean white noise level aps, 
computed in the range 0.5-1.5 mHz, reaches 3.06 x 10~^ 
m^s^^. With 1260 measurements, the velocity accuracy 
corresponds thus to 0.98 ms~^. 

Observed mode frequencies expected for solar-like 
stars correspond to low angular degree I and high radial 
order n and are well approximated by the asymptotic 
relation : 



with A:/ = iyn,l - Vn-l,l ■ 



(1) 



The large splitting /S.v reflects the average stellar 
density and e is a constant near unity (1.46 for the Sun) 
sensitive to the surface layers. In order to determine 
the average large splitting (Ai^), we calculate the comb 
respons e of th e power spectrum as proposed by Kjeldsen 
et al. ( 1995b ). An average comb response is computed 
using all peaks greater than 0.03 m^ s^^ (i.e. with 
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Table 2. Mode frequencies (in /iHz) of a Cen A estimated 
with the CLEAN procedure. The frequency resolution of 
our time series is 2.6 /iHz. 



Fig. 2. Power spectrum of the radial velocity measure- 
ments of a Cen A 




Fig. 3. Average comb response of the power spectrum of 
a Cen A computed with 0.5 /iHz resolution 



S/N greater than 3 in the amplitude spectrum) in the 
frequency range 1.7-3 mHz. The resulting comb response 
shown in Fig. 3 gives (Av) — 105.7 /xHz. 

In order to estimate the average amplitude of modes, 
we follow ed the m ethod described and used by Kjeldsen & 
Bedding ( |l995a|) . A simulated time series was generated 
with artificial signal plus noise using our window function. 
The generated signal corresponds to 20 modes centered at 





I = 


I = 1 


n=16 


1833.1 


- 


n=17 


- 


- 


n=18 


2041.1 


2095.6 


n=19 


2146.4 


2203.1 


n=20 


2251.0 


2308.4 


n=21 


2358.9 


2414.1 


n=22 


2463.3 


2519.1 


n=23 


2566.5 


2625.1 


n=24 


2671.6 


2732.7 


n=25 


- 


2838.0 


Aiy 


105.5 ± 0.6 


106.1 ± 0.4 



2.36 mHz, separated by Ai^/2 and convolved by a gaussian 
envelope. The generated noise is gaussian with a standard 
deviation of 1 ms~^. Several simulations show that the 
average amplitude of peaks is in the range 33-37 cms~^. 
An attempted identification of the strongest modes 
was processed with the CLEAN algorithm (Roberts et al. 
1987) in order to remove the effect of the window func- 
tion. The identified modes are listed in Table 2 and shown 
in Fig. 4. The values of n and I are deduced from the 
asymptotic relation (see Eq. m) assuming that the main 
peak corresponds to n = 21 and 1 = 0. With this config- 
uration, the parameter e is estimated to 1.35. This value 
is located in-between the Sun value (1.46) and the value 
deduced from theo retica l models (Guenther & Demarque 
^000| , Morel et al. ^000|) which give e in the range 1.07- 
1.27. The last row in Table 2 gives the average large split- 
ting for / = and 1 = 1. Both values are in agreement 
with this deduced from the comb response. 

4. Concluding remarks and prospects 

Our observations of a Cen A conduct to an obvious de- 
tection of p-modc oscillations. Several identifiable modes 
appear in the power spectrum between 1.7 and 3 mHz 
with an average large splitting of 105.7 /iHz and an am- 
plitude of about 35 cms^^. These characteristics are in 
full agreement with the expected values scaling from the 
Sun presented in section 1 and will bring constraints on 
theoretical models. 

This result demonstrates the power of the simultane- 
ous thorium RV method which reaches here the precision 
of 1 ms^^ at frequency higher than 0.5 mHz. We hope 
to obtain similar convincing seismological results on some 
bright solar-like stars with CORALIE. 

Observations of a Cen A were conducted by T. 
Bedding and his collaborators during the same period with 
the UCLES spectrograph on the 3.9-m Anglo Australian 
Telescope and the UVES spectrograph on the 8.2-m Very 
Large Telescope. We hope that these observations will con- 
duct to multi-site data combination in order to improve 
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Fig. 4. Identified p-mode oscillations of a Cen A in the power spectrum. Dashed lines and dotted lines correspond 
respectively to the modes I = and I — 1 listed in Table 2. 



the spectral window and explore further the P-mode spec- 
trum of a Cen A. 

We are aware that ground based observation is far from 
the expected accuracy of the future space missions like 
MOST, MONS and CORO T. Ho wever the future spectro- 
graph HARPS (Pepe et al. |200C| ), which will be installed 
on the 3.6-m ESO telescope at La Silla Observatory at 
the end of 2002, is expected to observe stars with magni- 



tude 5 dimmer than CORALIE (Bouchy et al. |2001|). It 
will thus be able to conduct complementary ground based 
asteroseismological study on a large sample of stars. 
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